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Abstract 

Given a polygon and a visibility range, the Myopic Watchman Problem with Discrete Vision (MW- 
PDV) asks for a closed path P and a set of scan points 5, such that (i) every point of the polygon is 
within visibility range of a scan point; and (ii) path length plus weighted sum of scan number along the 
tour is minimized. Alternatively, the bicriteria problem (ii') aims at minimizing both scan number and 
tour length. We consider both lawn mowing (in which tour and scan points may leave P) and milling 
(in which tour, scan points and visibility must stay within P) variants for the MWPDV; even for simple 
special cases, these problems are NP-hard. 

We show that this problem is NP-hard, even for the special cases of rectilinear polygons and Lea 
scan range 1, and negligible small travel cost or negligible travel cost. For rectilinear MWPDV milling in 
grid polygons we present a 2.5-approximation with unit scan range; this holds for the bicriteria version, 
thus for any linear combination of travel cost and scan cost. For grid polygons and circular unit scan 
range, we describe a bicriteria 4-approximation. These results serve as stepping stones for the general 
case of circular scans with scan radius r and arbitrary polygons of feature size a, for which we extend the 
underlying ideas to a 7r(^ + ^-j^) bicriteria approximation algorithm. Finally, we describe approximation 
schemes for MWPDV lawn mowing and milling of grid polygons, for fixed ratio between scan cost and 
travel cost. 

Keywords: Covering, Minimum Watchman Problem, limited visibility, lawn mowing, bicriteria 
problems, approximation algorithm, PTAS. 

1 Introduction 

Covering a given polygonal region by a small set of disks or squares is a problem with many applications. 
Another classical problem is finding a short tour that visits a number of objects. Both of these aspects have 
been studied separately, with generalizations motivated by natural constraints. 

In this paper, we study the combination of these problems, originally motivated by challenges from 
robotics, where accurate scanning requires a certain amount of time for each scan; obviously, this is also 
the case for other surveillance tasks that combine changes of venue with stationary scanning. The crucial 
constraints are (a) a limited visibility range, and (b) the requirement to stop when scanning the environment, 
i.e., with vision only at discrete points. These constraints give rise to the Myopic Watchman Problem with 
Discrete Vision (MWPDV), the subject of this paper. 

For a scan range that is not much bigger than the feature size of the polygon, the MWPDV combines two 
geometric problems that allow approximation schemes (minimum cover and TSP). This makes it tempting 
to assume that combining two approximation schemes will yield a polynomial-time approximation scheme 
(PTAS), e.g., by using a PTAS for minimum cover (Hochbaum and Maass [15|). then a PTAS for computing 
a tour on this solution. As can be seen from Figure [TJa) and (b), this is not the case; moreover, an optimal 
solution depends on the relative weights of tour length and scan cost. This turns the task into a bicriteria 
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problem; the example shows that there is no simultaneous PTAS for both aspects. As we will see in Sections U] 
and [5l a different approach allows a simultaneous constant- factor approximation for both scan number and 
tour length, and thus of the combined cost. We show in Section [71 a more involved integrated guillotine 
approach allows a PTAS for combined cost in the case of a fixed ratio between scan cost and travel cost. 




Figure 1: (a) An MWPDV solution with a minimum number of scans; (b) an MWPDV solution with a minimum 
tour length, (c) A minimum guard cover may involve scan points that are not from an obvious set of candidate 
points. 

A different kind of difficulty is highlighted in Figure [IJc), where an example called "wells" (12) is illus- 
trated: For a visibility range r that is large compared to the size of the niches (the "wells" ) , it may be quite 
hard to determine a guard cover of small size, since it is not clear what should be a small set of candidate 
guard locations that suffice for coverage. Each "well" is covered by a pair of points, whose locations need 
not correspond, e.g., to vertices in the arrangement of visibility polygons of the vertices of the polygon. In 
fact, because of this difficulty, there is no known nontrivial (with factor o{n)) approximation for minimum 
guard cover in simple n-gons. With an assumption about a sufficient set of guard points (e.g., vertices of the 
polygon or a grid within the polygon), an 0(log OPr)-approximation by Efrat and Har-Peled |10| is known. 
Also note that the optimal solution for our problem may change significantly with the relative weights be- 
tween tour length and scan cost: If the tour length dominates the number of scans in the objective function, 
an optimal tour can be forced to follow the row of niches on the right in the figure. We will show in Section 
6 how to obtain a constant- factor approximation for a bounded value -, a being the minimum side length 
of P. 

Related Work. Closely related to practical problems of searching with an autonomous robot is the 
classical theoretical problem of finding a shortest watchman tour, e.g., see [7j|8|. Planning an optimal set 
of scan points (with unlimited visibility) is the art gallery problem [20j . Finally, visiting all grid points of a 
given set is a special case of the classical Traveling Salesman Problem (TSP); see |16| . Two generalizations of 
the TSP are the so-called lawn mowing and milling problems: Given a cutter of a certain shape, e.g., an axis- 
aligned square, the milling problem asks for a shortest tour along which the (center of the) cutter moves, such 
that the entire region is covered and the cutter stays inside the region at all times. Clearly, this takes care of 
the constraint of limited visibility, but it fails to account for discrete visibility. At this point, the best known 
approximation method for milling is a 2.5-approximation [3]. Related results for the TSP with neighborhoods 
(TSPN) include I19j: further variations arise from considering online scenarios, either with limited vision 
[5] or with discrete vision [T31 [H] , but not both. The discrete visibility is intrinsic to the art gallery problem, 
but no tour is considered here. For this problem neither constant-factor approximation algorithms nor exact 
solution methods are known, recent results include an algorithm based on linear programming that provides 
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lower bounds on the necessary number of guards in every step and — in case of convergence and integrality — 
ends with an optimal solution by Baumgartner et al. [4]. Finally, [1] consider covering a set of points by a 
number of scans, and touring all scan points, with the objective function being a linear combination of scan 
cost and travel cost; however, the set to be scanned is discrete, and scan cost is a function of the scan radius, 
which may be small or large. 

For an online watchman problem with unrestricted but discrete vision, Fekete and Schmidt |14j present 
a comprehensive study of the milling problem, including a strategy with constant competitive ratio for 
polygons of bounded feature size and with the assumption that each edge of the polygon is fully visible from 
some scan point. For limited visibility range, Wagner et al. |23j discuss an online strategy that chooses an 
arbitrarily uncovered point on the boundary of the visibility circle and backtracks if no such point exists. 
For the cost they only consider the length of the path used between the scan points, scanning causes no cost. 
Then, they can give an upper bound on the cost as a ratio of total area to cover and squared radius. 

Our Results. On the positive side, we give a 2.5-approximation for the case of grid polygons and a 
rectangular range of unit-range visibility, generalizing the 2.5-approximation by Arkin, Fekete, and Mitchell 
[3] for continuous milling. The underlying ideas form the basis for more general results: For circular scans 
of radius r = 1 and grid polygons we give a 4-approximation. Moreover, for circular scans of radius r 
and arbitrary polygons of feature size a, we extend the underlying ideas to a. Tr{^ + ^^-ji )-approximation 
algorithm. All these results also hold for the bicriteria versions, for which both scan cost and travel cost 
have to be approximated simultaneously. Finally, we present a PTAS for MWPDV lawn mowing and a 
PTAS for MWPDV milling, both for the case of fixed ratio between scan cost and travel cost. 

The rest of the paper is organized as follows. In the following Section [2] we give the notation and 
formally define the Myopic Watchman Problem with Discrete Vision. Section [3] provides a NP-hardness 
proof. Approximation algorithms for grid polygons and rectangular unit scan range, grid polygons and 
circular unit scan range as well as for general polygons and circular scan range are presented in Sections |4l 
[5] and [6l respectively. A description of polynomial-time approximation schemes for both the lawn mowing 
and the milling variant are given in Section [T] In the final Section |8] we discuss possible implications and 
extensions. 

2 Notation and Preliminaries 

We are given a polygon P. In general, P may be a polygon with holes; in Sections 3, 4 and 5, P is an 
axis-parallel polygon with integer coordinates. 

Our robot, R, has discrete vision, i.e., it can perceive its environment when it stops at a point and 
performs a scan, which takes c time units. From a scan point p, only a ball of radius r is visible to R, either 
in Loa- or L2-metric. A set S of scan points covers the polygon P, if and only if for each point q Cz P there 
exists a scan point p G S such that q sees p (i.e., qp C P) and \qp\ < r. 

We then define the Myopic Watchman Problem with Discrete Vision (MWPDV) as follows: Our goal is 
to find a tour T and a set of scan points S{T) that covers P, such that the total travel and scan time is 
optimal, i.e., we minimize t(T) = c ■ \S{T)\ + L(T), where L{T) is the length of tour T. Alternatively, we 
may consider the bicriteria problem, and aim for a simultaneous approximation of both scan number and 
tour length. 

3 NP-Hardness 

Even the simplest and extreme variants of MWPDV lawn mowing arc still generalizations of NP-hard prob- 
lems. 

Theorem 1. 

(1) The MWPDV is NP-hard, even for polyominoes and small or no scan cost, i.e., c <^ 1 or c = 0. 

(2) The MWPDV is NP-hard, even for polyominoes and small travel cost, i.e., c 3> 1. 
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Figure 2: Top row: The polygonal piece for a variable and the connection of an edge corridor. Bottom row: a 
placement corresponding to "true" (left), and a placement corresponding to "false" (right). The light gray scan is 
used in both cases. 



(3) The MWPDV is NP-hard, even for polyominoes and no travel cost, i.e., t(T) = \S\. 

Proof. The first claim is a result of the hardness of minimum cost milling, see [5]. The second claim is an 
easy consequence of the NP-hardness of Hamiltonicity of Grid Graphs (HGG) [16]: Given an instance 
G of HGG with n vertices, turn it into an instance of MWPDV by scaling the grid graph G by a factor of 
two, and replacing each grid point of G by a 2x2-square. This yields a canonical set of n scan points that 
is contained in any optimal WMPDV tour; traveling these with a tour length 2n is possible if and only the 
graph G is Hamiltonian. 

The third claim is closely related to a minimum cover problem by visibility discs; however, the MWPDV 
requires that the scan points must be inside of the polygonal region. We give a proof along the lines of [5] , 
based on a reduction of the NP-hard problem Planar 3SAT, a special case of 3SAT in which the variable- 
clause incidence graph H is planar. As a first step, we construct an appropriate planar layout of the graph 
H, e.g., by using the method of Rosenstiehl and Tarjan |21j . This layout is turned into a grid polygon by 
representing the variables, the clauses and the edges of G. An example for the variable component is given 
in Figure m 

The variable gadgets allow two ways for locating a minimum number of scan points. The first (the black 
points in Figure [2]) relates to a setting of "true", the other (the circles in Figure [2]) to a setting of "false". For 
the variable setting "true" the scan squares are pushed further into the edge corridor. These edge corridors 
arc similar to the variable-circles — bendings are done accordingly. In order to have the same number of 
points and circles, edge corridors may be added, that do not end in another polygonal piece, but assure this 
parity (with circles at the edge corridor). 

A clause component is given in Figure [H Edge corridors of the three associated variables (each with the 
appropriate truth setting) meet in the polygonal piece for the clause (dark gray in Figure |3]). If and only 
if the clause is satisfied, i.e., if at least in one edge corridor a scan square is placed at a black point, three 
additional scans suffice to cover this polygonal piece. Otherwise, four scans are necessary. 
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Given the components defined above we can compute the parameter k, the number of scan squares 
necessary to cover the entire resulting polygon P. k is polynomial in the number of vertices of G and part of 
the input. All vertices of the resulting P have integer coordinates of small size, their number is polynomial 
in the number of vertices of G. This shows that the problem is NP-hard. □ 




(a) (b) (c) 



Figure 3: A clause component with the polygonal piece (dark gray) and the three according edge corridors (blue) (a). 
A placement corresponding to "true" (b), and a placement corresponding to "false" (c). The green squares are the 
scans necessary inside the clause. 

□ 

4 Approximating Rectilinear MWPDV Milling for Rectangular 
Visibility Range 

As a first step (and a warmup for more general cases), we give an approximation algorithm for rectilinear 
visibility range in rectilinear grid polygons. 

The following lemma allows us to focus on visiting and scanning at grid points. 

Lemma 2. For a rectilinear pixel polygon P there exists an optimum myopic watchman tour T* such that 
all scan points are located on grid points: 

3T*V(x„ y,) e S{T*) : x, e Z A £ Z. (1) 

Proof. Let T be an optimal tour, with scan points not located on grid points. Consider the vertical and 
horizontal strips of pixels in P of maximal length. W.l.o.g., we start with the horizontal strips. For every 
strip, we shift the scans, such that the x-coordinates are integers (starting from the boundary, i.e., with 
distance 1 to the boundary if possible, and away from non- reflex corners). The tour will not be longer, we 
cover not less; in case we are able to reduce the number of scans per strip by one we have a contradiction to 
t being optimal. After applying this to all horizontal strips we proceed analogously for the vertical strips. 
Hence, we have an optimal tour, with all scan points located on grid points. □ 

Our approximation proceeds in two steps: 

(I) Construct a set of scan points that is not larger than 2.5 times a minimum cardinality scan set. 

(II) Construct a tour that contains all constructed scan points and that does not exceed 2.5 times the cost 
of an optimum milling tour. 

The main idea for the second step is based on the 0{n log n)-time 2.5-approximation algorithm for milling 
from Arkin et al. [3]. The resulting tour consists of three parts, see Figure 2] for an example, Lqpt being 
the optimal milling tour length: 
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(1) a "boundary" part: B C P is the inward offset region of all points within P that are feasible placements 
for the center of the milling cutter. For a milling problem, B is connected. Tracing the boundary SB 
of B, let PsB denote the region milled by this route. {SB may not be connected (if P features holes), 
the pieces arc SBi.) The length of SB, Lsb is a lower bound on Lqpt- 

(2) a "strip" part: Pint P\Psb — if nonempty — can be covered by a set of k horizontal strips E,;. The 
y-coordinates of two strips differ by multiples of 2. Then, let L^tr = X)i=i and this is again a lower 
bound on the length of an optimal milling tour: Lqpt ^ Lgtr- 

(3) a "matching" part: the strips and the boundary tour have to be combined for a tour. For that purpose, 
consider the cndpoints of strips on SBf. every SBi contains an even number of such cndpoints. Hence, 
every SBi is partitioned into two disjoint portions, Mi{SBi) and M2{SBi). Using the shorter of these 
two {A'h{SBi)) for every SBi we obtain for the combined length, Lm'- Lm < Lgtr/'^ < Lopt/'^- 

The graph with cndpoints of strip lines plus the points where strip line touches a SBi as vertices is connected 
by three types of edges — the center lines of the strips, the SBiS and the M^,{SBi)s. Every vertex has degree 
4, hence, an Eulerian tour gives a feasible solution. 
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Figure 4: A polygon P, consisting of pixels. The parts for the approximative milling tour are indicated: region Psb., 
is shaded (light blue), SBi is highlighted in blue. 

Now we describe how to construct a covering set of scan points: 
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1. Let 5*46 be the "even quadruple" centers of all 2x2-squares that are fully contained in P, and which 
have two even coordinates. 

2. Remove all 2x2-squares corresponding to 5*46 from P; in the remaining polyomino Pie, greedily pick a 
maximum disjoint set S^o of "odd quadruple" 2x2-squares. 

3. Remove all 2x2-squares corresponding to S^o from P^e] greedily pick a maximum disjoint set S3 of 
"triple" 2x2-squares that cover 3 pixels each in the remaining polyomino P4e.4o, 

4. Remove all 2x2-squares corresponding to 5*3 from P4e,4o; hi the remaining set P4e,4o,3 of pixels, no three 
can be covered by the same scan. Considering edges between pixels that can be covered by the same 
scan, pick a minimum set of ("double" S2 and "single" 5"!) scans by computing a maximum matching. 




Figure 5: An example for our approximation method: The set of "even quadruple" scans is shown in grey; the "odd 
quadruple" scans are green. A possible (greedy!) set of "triple" scan is shown in blue, leaving the maximum matching 
(and the corresponding "double scans") shown in red. The leftover single pixels are yellow. The ellipse indicates a 
part that is covered by three scans instead of two; the triple scan with adjacent single and double scans could be 
covered by two triple scans. 



Claim 3. The total number of scans is at most 2.5 times the size of a minimum cardinality scan set. 

Claim 4. All scan points lie on a 2. 5 -approximative milling tour. 

For the first claim, let Smin be the size of a cardinality scan set. Observe that 54e U 6*40 corresponds to a 
set of disjoint 2x2-squares that are fully contained in P; hence, it follows from a simple area argument that 

\Sie U 5*40! ^ Smin- 

When considering the set of pixels in /'4e.4o, not more than three can be covered by the same scan; in 
P4e,4o,3 we compute an optimal solution. This implies that the only way to get a smaller cover in P4e.4o is 
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to change the choice of triple scans; this means that the pixels covered by a triple scan have to be allocated 
differently. Taking into account that Pie,io does not contain any 2x2-squares, a simple case analysis shows 
that only two possible improvements are possible: 

• replacing a triple, a double and a single by two triples (as shown in the left part of Figure [5]) , or 

• replacing a triple and two singles by a triple and a double. 

With s^jj^ being the minimum number of scans required for covering P4e,4o, we get 1 5*3 U 5*2 U | < fsf^jn < 
fsmin- In total wc get a solution with not more than |smin scans. 

For the second claim, we use a milling tour constructed as in [3] and described above: Choosing the 
strips to be centered on even y-coordinates allows us to visit all scan points in S^e- Clearly, all pixels in P^^ 
are adjacent to the boundary of other scans involve boundary pixels, allowing them to be visited along the 
"boundary" part. (One minor technical detail is shown in Figure |6j In order to visit the center of a triple 
scan, we need to reroute the boundary part of the tour to run through a reflex vertex; this does not change 
the tour length.) 



Figure 6: Modifying the boundary part of W{t): at a reflex vertex, the gray path is used instead of SBi. 

This concludes the proof. We summarize: 

Theorem 5. A polyomino P allows a MWPDV with rectangular vision solution that contains at most 2.5 
times the minimum number of scans necessary to scan the polygon, and has tour length at most 2.5 times 
the length of an optimum milling tour. 



5 Approximating Rectilinear MWPDV Milling for Circular Visi- 
bility Range 

In this section we give an approximation algorithm for MWPDV milling in case of a circular range of visibility. 
Again, P is a polyomino, we consider r = 1 and the tour length is measured according to the Li-distance. 
The approximation is for the case of no given starting point or in case of a given starting point located on a 
grid point (in the polyomino or on its boundary). 

When considering a circular scan range, one additional difficulty arc boundary effects of discrete scan 
points: While continuous vision allows simply sweeping a corridor of width 2r by walking down its center 
line, additional cleanup is required for the gaps left by discrete vision; this requires additional mathematical 
arguments, see Figure [T] 

We overlay the polyomino with a point grid as in Figure [3 left, i.e., a diagonal point grid with L2-distance 
of in-between points, and use all points of this grid that coincide with a grid point in P (Figure O left). 
For the tour T our exploration strategy starts at a boundary grid point, proceeding counterclockwise along 
the boundary and taking a scan at every point of the overlayed grid located on the boundary. Only using 
this we would end up with a tour covering an area around the boundary. 
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Figure 7: A finite number of scan points along tlie center line of a corridor of width 2r does not suffice to cover the 
corridor with a circular scan range of radius r. 
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Figure 8: Left: Point grid (light gray) with grid points within a polyomino (black) in dark gray. Right: Circular 
visibility ranges of the grid points covering the plane, one square of side length \/2 is indicated in yellow. 



For the movement in-between interior scan points and the boundary we use horizontal strips located on 
grid lines (and distance 1 to the boundary), cp. Figure IH] for an example. In order to combine these grid lines 
with the boundary path for a tour, we link strips to the left boundary. Two of these will in general be linked 
on the right-hand side. In case there is an odd number of strips between the upper and lower boundary of 
P, the scan points located on the bottommost strip are visited by using a path of (Li-)length 2, from the 
boundary or strips with another y-coordinate of the leftmost point, cp. Figure|9l down right. The strips that 
get linked are always determined by the leftmost boundary. In case there are other parts of P that have a 
left boundary (that is vertical edges with polygon to the right and the exterior to the left) whose cardinality 
of strip lines differ by an odd number, scan points from the topmost strip arc linked by two vertical steps of 
length 1 to the upper boundary, cp. Figure [9l top right. 

Hence, we yield a closed tour (linking always two strips we always end up at the left boundary, the rest 
is a tour along the boundary with small loops of length 2). We still need to show that P is covered and have 
to consider the competitive ratio of our strategy. Let T be the tour determined by our strategy, and T* be 
an optimal tour. 

Claim 6. The scan points the strategy positions cover P. 

Proof. Including the scan points on the boundary the entire polyomino is covered by the scan point grid 
defined above: The scan range being a i2-circle, these scans cover the interior of the polyomino — each circle 
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Figure 9: A polyomino P with the tour given by our strategy. Scan points are displayed in black. The horizontal 
strips of total length Lstrips are indicated in yellow, the tour is in blue for the links to the boundary and in-between 
strips (solid) as well as for connections of points (dash-dotted), parts located on the strips are indicated in red. 
(In-between those parts the tour runs on the boundary.) 

covers at least a square of side length \/2 (with the vertices located on grid points that do not belong to the 
overlayed grid), cp. Figure [3 right. □ □ 

Lemma 7. L(T) < 4 • L{T*) + 8 

Proof. Our tour consists of the tour along the boundary of P, the strips and connections. Let ibounrf be the 
length of the boundary, i.e., P's perimeter. Moreover, we define Si: 

• to be Ssi (cp. Section |4l) whenever the corresponding cutter fits into P, 

• in-between these Sbi a shortest path on the grid, 

• plus (for corridors of P of width 1) shortest path from the SbiS such that every pixel of P is visited by 
a part of 5i (a pixel is visited if one of its vertices or edges is reached). 

Finally, let the length of all strips be Lstrips ■ 

We need to show that these three elements cover the length of our chosen path. Movements along the 
boundary are covered by Lbound, movements along strips by Lstrips- The right-hand side links of two adjacent 
strips lie on 6i, the left-hand side connections can be turned by 90° (only for charging) and lie on Si as well. 
The scan point connecting paths of length 2 can be swung open and lie on the strip of the relative (not 
used) strip or one part (length 1) on 6i for interior strip points or endpoints of strips, respectively. Thus, 
the combined lengths of these are an upper bound for our tour length. 

Using the fact that for an orthogonal polygon of n vertices, r of which are reflex, n = 2r + 4 holds, see 
[20] . Lemma 2.12, we have Ltound = Lsi + 8. Furthermore, with L^tr as defined as in Section |4l we have 
Lstrips < 2 • Lstr < 2 • L{T*) . Finally, the milling argument shows L^i < L{T*). Altogether, we have: 

L(T) < Ltound + Lstrips + Lsi < 2 • Lsi + 8 + Lstrips < 4 • L{T*) +8. □ □ 
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Lemma 8. \S{T)\ < 4 • \S{T*)\ 

Proof. Let N{P) be the number of pixels of a polyomino P. Moreover, let V{N{P)) be the maximum ratio 
||^^ljrij| for all polyominos with N{P) pixels. Kershner [T7](cp. [22]) showed (with the notation introduced 

above) that D{r) = "'"'^jfiJP^ > that is, for r = 1: |5(r*)| > ^ • N{P). The placement of scan 

points on the diagonal grid allows us to bound the number of scan points by 2 + {N{P) — 1) = N{P) + 1- If we 
define F{N{P)) :— ^^^"^^^ , this gives us an upper bound on V{N{P)). F{n) is monotonically decreasing 

in n. Furthermore, Fin) = 4 ^ n = _3v^ « 1.85322, that is: 

' ^ ' 8-3V3 

V{N{P)) < F{N{P)) < 4 for N{P) > 2 (2) 

(For N{P) = 1 the optimum needs at least one scan, we need at most 2.) □ □ 

Theorem 9. A polyomino P allows a MWPDV solution for a circular visibility range with r = 1 that is 
A-competitive. 

Proof. With Lemmas [7] and [8] we have: 



t{T) = c-\SiT)\+L{T) 

< c-4- |5(T*)| +4-i(T*) + 
= 4-(c-|5(T*)|+L(r*)) + i 



□ □ 



6 Approximating General MWPDV Milling for a Circular Visi- 
bility Range 

In this section we discuss MWPDV milling for a circular visibility range r in general polygons. As discussed 
in Section 1, even the problem of minimum guard coverage has no known constant-factor approximation; 
therefore, we consider a bounded ratio r/a between visibility range and feature size, i.e., minimum side 
length. 

Just as in the rectilinear case for a rectilinear scan range, see Section |4l our approximation proceeds in 
two steps: 

(I) Construct a set of scan points that is within a constant factor of a covering set of minimum cardinality. 

(II) Construct a tour that contains all constructed scan points and is within a constant factor of the cost 
of an optimum milling tour. 

We start with a description of the second step, which will form the basis for the placement of scan points. 
Just as in the rectilinear case, we consider three parts. 

(1) A "boundary" part: Above we described tracing 6B the boundary of B, causing a tour length of Lsb- 
Here, we use two "boundary tours" within distance of (at most) ^r and (at most) ^r to the boundary, 
TRl and TR2 of length Ltri and Ltr2 , respectively. Then, we have: 

Ltbi + Ltr2 = 2-Lsb<2-L{T*) (3) 

(The length of the three tours differs at the vertices: drawing a line perpendicular there from TR2 to 
TRl the Intercept Theorem shows that the distance to the diagonal through the vertices of all tours 
on TRl is twice as much as on the boundary tour with distance r to the boundary.) 

The two "boundary" tours allow us to cover a corridor of width 2r with a bounded number of scans, 
while ([3]) enables us to bound the tour length in terms of the optimal length. 
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(2) A "strip" part: For the interior we use strips again: Pint '■= P\Psb — if nonempty — can be covered 
by a set of fci horizontal strips Y,}. The y-coordinates of two strips differ by muhiples of 2r. We can 
consider another set of strips, E^, shifted by r. Then, let L^^^ = Xli^i -^s^.- Similar to the argument 
for Loo, we have 1]^^ + <2- L{T*). 

(3) A "matching" part: In order to combine the two "boundary parts" and the two sets of strips for a tour 
we add two more set of sections. 

— The center lines of the strips have a distance of r to the boundary, thus they do not yet touch 
TRl. Consequently, we add l/2r to each center line (on each end). For that purpose, we consider 
the matchings as defined above. (Consider the endpoints of strips on 6Bi: every 5Bi contains 
an even number of such endpoints. Hence, every 5Bi is partitioned into two disjoint portions, 
Mi{5Bi) and M2{SBi). Using the shorter of these two {M^.{6Bi)) for every 5Bi we obtain for the 
combined length, Lm- Lm < Lstr/2 < L(T*)/2.) Because two strips are at least a distance of r 
apart, the connection to TRl costs less than 1/2 • Lm < 1/2 • Lstr/2 < L{T*)/A. 

— Moreover, we consider the above matchings defined on TRl and insert the shorter sections of 
the disjoint parts, (MliSBi)), for every 5Bi. The Intercept Theorem in combination with the 
analogously defined sections on TR2 enables us to give an upper bound of L^/i < Lstr < L{T*). 

Starting on some point on TRl, tracing the strips, and the inner "boundary" TR2 at once when passing 
it yields a tour; the above inequalities show that L{T) < 21/4 • L{T*). 

Now we only have to take care of (I), i.e., construct an appropriate set of scan points. For the "boundary" 
part we place scans with the center points located on TRl and TR2 in distance -y/S • r (see Figure [T0|) if 
possible, but at corners we need to place scans, so the minimum width we are able to cover with the two 
scans (on both tours) is a. For the "strip" part the distance of scans is also • r on both strip sets, exactly 
the distance enabling us to cover a width of r, see Figure [TOl 



Figure 10: The dashed lines indicate a length of r, so the dash-dotted line has a length of -\/3 ■ r, the distance of two 
scans. 

It remains to consider the costs for the scans. We start with the inner part. Taking scans within a 
distance of -y/S • r, we may need the length divided by this value, plus one scan. We only charge the first part 
to the strips, the (possible) additional scans arc charged to the "boundary" part, as we have no minimum 
length of the strips. The optimum cannot cover more than irr^ with one scan. Let Lstr = '^''^^{Ll^j,, L^tr)' 

\siT )| > \S{T)\ < ^ ^ < ^ - ^ (4) 

Finally, we consider the "boundary". We assume Lsb > 1. So |iS(T*)| > We may need to scan within 

a distance of a — on two strips — , need additional scans and have to charge the scans from the "strip" part, 
hence, this yields: \S{T)\ < + 1 + Consequently, for r > a: 

\S{T)\ ^ ^ ^ _^ ^ (.5^ 



\S{T*)\ - a 2 2 
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Theorem 10. A polygon P allows a MWPDV solution that contains at most a cost of max{'^ , + + f) 
times the cost of an optimum MWPDV solution (for r > a). 

Note that Theorem [TOl covers the case from Section [S] however, instead of the custom-buih factor of 4 it 
yields a factor of 2 • tt. 

7 A PTAS for MWPDV Lawn Mowing 

We describe in detail here the following special case, and then discuss how the method generalizes. Consider 
a polyomino P (the "grass") that is to be "mowed" by a fc x fc square, M. At certain discrete set S{T) of 
positions of M along a tour T, the mower is activated (a "scan" is taken), causing all of the grass of P that 
lies below M at such a position to be mowed. For complete coverage, we require that P be contained in the 
union of k x k squares centered at points S{T). Between scan positions, the mower moves along the tour T. 

In this "lawn mower" variant of the problem, the mower is not required to be fully inside P; the mower 
may extend outside P and move through the exterior of P, e.g., in order to reach different connected 
components of P. Since P may consist of singleton pixels, substantially separated, the problem is NP-hard 
even for k = 1, from TSP. 

Here we describe a PTAS for the problem. We apply the m-guillotine method, with special care to handle 
the fact that we must have full coverage of P. Since the problem is closely related to the TSPN |9l [19] , we 
must address some of the similar difficulties in applying PTAS methods for the TSP: in particular, a mower 
centered on one side of a cut may be responsible to cover portions of P on the opposite side of the cut. 

At the core of the method is a structure theorem, which shows that we can transform an arbitrary tour 
T, together with a set S{T) of scan points, into a tour and scan-point set, {Tq, S{Ta)), that are m-guillotine 
in the following sense: the bounding box of the set oi k x k squares centered at S{T) can be recursively 
partitioned into a rectangular subdivision by "m-perfect cuts". An axis-parallel cut line i is m-perfect if 
its intersection with the tour has at most m connected components and its intersection with the union of 
k X k disks centered at scan points consists of m disks or "chains of disks" (meaning a set of disks whose 
centers lie equally spaced, at distance fc, along a vertical/horizontal line); see Figure [TT] (The definition of 
m-perfect in |18[ has a slightly different specification, in terms of the number of endpoints of the connected 
components of Tn ^, but it is within a constant factor equivalent to what we define here.) 

The structure theorem is proved by showing the following lemma: 

Lemma 11. For any fixed m = [l/el and any choice of (T,S{T)), one can add a set o/ doubled bridge 
segments, of total length 0{\T\/m), to T and a set of 0{\S{T)\/m) bridging scans to S{T) such that the 
resulting set, {Tg,S{Tq)), is m-guillotine, with points S (Tg) on tourTo and withTc containing an Eulerian 
tour of S{Tg). 

Proof. For a vertical line, Ix, through coordinate x, let /i(x) denote the length of the m-span of Ix with 
respect to R and T: fi{x) = if, within i?, tx intersects T in at most 2m connected components; otherwise, 
if Ix intersects T H i? in components ci,C2, . . . ,Cif, K > 2m, then fi{x) is the distance (along ix) from 
component Cm to component CK-m+i- Similarly, we define the length, f2{x), of the m-scan-span of ix with 
respect to R and S{T): f2{x) = if, within R, ix intersects at most 2m of the scan disks {k x k squares) 
centered at points S{T); otherwise, if ix intersects scan-disks Di,D2, ■ ■ ■ ,Dk, K > 2m, then f2{x) is the 
distance (along ix) from disk D„i to disk DK-m+i- We think of fi{x) as the "cost" to construct a vertical 
bridge for T at position x, and f2{x) as the "cost" to add a sequence (chain) of scan disks, centered along 
£x, and a detour of the tour (of length 0{k) per disk) that visits their centers. (We similarly define costs 
gi{y) and 52(2/) for constructing bridges along a horizontal cut iy through coordinate y.) The total cost 
associated with a vertical cut at £x is proportional, then, to /i(x) -I- /2(a;), and the cost of a horizontal cut 
at £y is proportional to giiy) + 32(2/)- (Here is where we are using the fact that the total cost is a linear 
combination of tour length and number of scans, with a fixed bound, c, on the relative cost ratio of length 
versus number of scans and that the scan disk has constant size k.) 

In order to charge off the cost of constructing bridges along the m-span and adding a chain of scan disks 
(and a subtour linking them) along the m-scan-span, we use the notion of "chargeable length" based on the 
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Figure 11: The green axis-parallel cut line t intersects the tour T (shown in red) in 4 connected components and 
intersects the scan ranges (shown in blue), centered at scan points, in a single chain of fc x fc disks. (Here, k = 2.) 

"m-dark" length of ix H i?, with a notion of "m-darkness" determined not just from the tour T, but also 
from the set of scans S{T). Specifically, a subset ab oi l^OR is said to be m-dark with respect to T if for any 
p e ab, the rightwards and leftwards rays from p each cross at least m (vertical) segments of T before exiting 
R. Similarly, ab is said to be m-dark with respect to the scan disks S(T) if for any p e ab, the rightwards 
and leftwards rays from p each intersect at least m scan disks centered on S{T) before exiting R. If a cut 
is made along £x, then the m-dark with respect to T portion of the cut can be charged off to the left /right 
sides of segments of T lying to the right/left of £x, distributing the charge to be (l/m)th to each of the m 
segments first hit. Similarly, the portion of the cut that is m-dark with respect to S{T) can be charged off 
to the scan-disks (or, more precisely, to their total perimeter, which is proportional (via constant k) to their 
cardinality) . 

The key observation, then, is that there must exist a "favorable" vertical cut or horizontal cut £y 
for R such that the chargeable length of the cut is at least as long as the cost of the cut. This follows 
from the usual argument (|18]). using the fact that Jx^uifii^) + f2{x))dx = ^^^j^h{y)dy, where h{y) is 
the chargeable length associated with the horizontal cutt ly, and assuming, without loss of generality, that 
Jx(zR{fi{x) + f2{x))dx > Jy^j^igiiy) + g2{y))dy: There must exist a value y* where gi{y*) + g2{y*) < h{y*), 
which then defines a favorable cut iy- for which the cost of constructing the m-span bridge and the m-scan- 
span sequence of scan disks is chargeable to lengths of T and disks of S{T) in such a way that no length of 
disk gets charged more than an amount proportional to (l/m)th of its length/count. 

Once a favorable cut is found with respect to one rectangle R, the cut partitions the problem into two 
subrectangles, and the argument is recursively applied to each. The end result is an m-guillotine subdivision 
of the original bounding box of T. □ □ 
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Figure 12: An example for the boundary information witliin a rectangle R. The "portals" are depicted in violet, the 
bridges in green, the disk bridges in yellow, the scan positions such that a k x k square centered at each position 
intersects the corresponding side of R in brown and the connection pattern in turquoise. (The k x k scan range is 
shown by dotted lines.) 

The algorithm is based on dynamic programming to compute an optimal m-guillotine network. A sub- 
problem is specified by a rectangle, R, with integer coordinates. The subproblem includes specification of 
boundary information, for each of the four sides of R; see Figure [T^ The boundary information includes: (i) 
0{m) integral points ("portals") where the tour is to cross/touch the boundary, (ii) at most one (doubled) 
bridge and one disk-bridge (chain) per side of R, with each bridge having a parity (even or odd) specify- 
ing the parity of the number of connections to the bridge from within R, (iii) 0{m) scan positions (from 
S{T)) such that a. k x k square centered at each position intersects the corresponding side of R, (iv) a 
connection pattern, specifying which subsets of the portals/bridges are required to be connected within R. 
There are a polynomial number of subproblems. For a given subproblem, the dynamic program optimizes 
over all (polynomial number of) possible cuts (horizontal or vertical), and choices of bridge, disk-bridge, 
parity assignments, and compatible connection patterns for each side of the cut. The result is an optimal 
m-guillotine network, with doubled bridges (so that it contains an Eulerian subgraph spanning the nodes), 
and a scan-point set, S, visited by the network, such that the union oi k x k squares centered at the points 
S covers the input polygon P. Since we know, from the structure theorem, that an optimal tour T, together 
with a set S{T) of covering scan points can be converted into a tour and scan-point set, {Tq,S{Tq)), that 
are m-guillotine, and we have computed an optimal such structure, we know that the tour we extract from 
our computed network approximates optimal. We summarize: 

Theorem 12. There is a PTAS for MWPDV lawn mowing of a (not necessarily connected) set of pixels by 
a k X k square. 

The Milling Variant. Our method applies also to the "milling" variant of the MWPDV, in which the 
scans all must stay within the region P, provided that P is simple (no holes), and all of P is reachable by 
a scanner (fc x k square). We describe the changes that are required. Subproblems are defined, as before, 
by axis-aligned rectangles R. The difficulty now is that the restriction of i? to P means that there may be 
many {^{n)) vertical/horizontal chords of P along one side of R. We can ignore the boundary of P and 
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construct an m-bridge (which we can "afford" to construct and charge off, by the same arguments as above) 
for T, but only the portions of such a bridge that he inside P (and form chords of P) are traversable by our 
watchman (since the other portions are outside P). For each such chord, the subproblcm must "know" if the 
chord is crossed by some edge of the tour, so that connections made inside i? to a chord are not just made to 
a "dangling" component. We cannot afford to specify one bit per chord, as this would be 2^*^"^ information. 
However, in the case of a simple polygon P, no extra information must be specified for the subproblcm - a 
chord is crossed by T if and only if the mower (scan) fits entirely inside the simple subpolygon on each side of 
the chord; if the subpolygon outside of R, on the other side of a chord, does not contain a kx k square, then 
we know that the entire subpolygon is covered using scanned centered within R. Exploiting this fact, the 
dynamic programming algorithm of our PTAS is readily modified to the MWPDV milling problem within a 
simple rectilinear polygon. 

Theorem 13. There is a PTAS for MWPDV milling of a simple rectilinear polygon by a k x k square. 

8 Conclusion 

A number of open problems remain. Is it possible to remove the dependence on the ratio (r/a) of the 
approximation factor in our algorithm for general MWPDV milling? This would require a breakthrough for 
approximating minimum guard cover; a first step may be to achieve an approximation factor that depends 
on log(r/a) instead of {r/a). 

For combined cost, we gave a PTAS for a lawn mowing variant, based on guillotine subdivisions. The 
PTAS extends to the milling case for simple rectilinear polygons. Wc expect that the PTAS extends to 
other cases too (circular scan disks, Euclidean tour lengths), but the generalization to arbitrary domains 
with (many) holes seems particularly challenging. Our method makes use of a fixed ratio between scan cost 
and travel cost; as discussed in Figure 1, there is no PTAS for the bicritcria version. 
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